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Abstract 

Let / be a positive smooth function on a close Riemann surface (M,g). The / — energy 
of a map u from M to a Riemannian manifold (iV, h) is defined as 



E f (u) = [ f\Vu\ 2 dV a 



In this paper, we will study the blow-up properties of Palais-Smale sequences for Ef. We will 
show that, if a Palais-Smale sequence is not compact, then it must blows up at some critical 
points of /. As a sequence, if an inhomogeneous Landau-Lifshitz system, i.e. a solution of 

u t — u x Tf(u) + Tf(u), u : M — > S 2 

blows up at time oo, then the blow-up points must be the critical points of /. 
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1 Introduction 

Let (M, g) and (N, h) be two Riemannian manifolds. A C 1 -smooth map u from M into N is 
called a harmonic map if and only if u is a critical point of the energy functional E(v), which 
is defined in local coordinates by 



E{v) = [ Trace g (v* h)dV g , 

J M 



where 

■ ■ du a du^ 
Trace g (v*h) = g»——h ol p{u). 

It is well-known that the energy functional is conformally invariant when dim(M) = 2. 

In this paper we would like to study a class of C 1 -smooth maps from a Riemann surface into 
a compact Riemannian manifold which are defined as the critical points of the inhomogeneous 
energy functional written as 



E f (v)= I Trace g (v*h)fdV g , 

J M 
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where / is a smooth real function. In [L] and |E-L| (see page 48, (10.20)), such maps are 
called /-harmonic from M into N. Obviously, they are just harmonic maps if / = 1. Moreover, 

2 

when m = dim(M) ^ 2, an /-harmonic map is nothing but a harmonic map from (M, f n ~ 2 g) to 
(N, h). In local coordinates, the /-harmonic map satisfies the following Euler-Lagrange equation 

fr{u) + V/ • Vn = 0. 

Here r(u) is the tension field of n which can be written as 

To see the physical motivation for the /-harmonic maps, we consider a smooth domain Q in 
the Euclidean space M. m . An inhomogeneous Heisenberg spin system is given by 

d t u = f(u A An) + V/ • (u A Vn), 

where / is a real-valued function defined on Q, u(x,t) £ S 2 , A denotes the cross products in R 3 
and A is the Laplace operator on M. m . Physically, the function / is called the coupling function, 
and is the continuum limit of the coupling constants between the neighboring spins. It is easy to 
see that if u is a smooth stationary solution of the above equation, then u is just an /-harmonic 
map from Q into S 2 . Indeed, in this case the tension field of u can be written as An + |Vn| 2 n, 
therefore, the right hand side of the above equation can be expressed by u A (/r(u) + V/ • Vn), 
and u satisfies the following equation 

fr(u) + V/ • Vn = 0. 

The above inhomogeneous Heisenberg spin system is also called inhomogeneous Landau- 
Lifshitz system. Landau and Lifshitz also suggested considering the following dispersive system 

d t u = u A (/r(n) + V/ • Vn) - u A (u A (/r(u) + V/ • Vn)), 

with an initial value condition 

n(0) = n . 

For the well-known equation, Tang [Tj proved that it admits a global weak solution which is 
smooth except for finitely many points, if the domain manifold M is 2-dimensional closed , / is 
a smooth positive function and the initial value map belongs to W 1 ' 2 (M,S 2 ) (see also fSfr and 
|G-H| ) . The bubbles which the weak solution blows are called as the magnetic bubbles ( |Shj ) . 
A natural question arises, where do the bubbling points of the Landau- Lifshitz equation locate? 
In this paper, we intend to answer this problem partially. 

Throughout this paper, we will always assume that / is smooth and positive. In order to 
answer the above question, mathematically we need to consider the convergence and bubbling 
of the sequence of /-harmonic maps with coupling function /. Precisely, we obtain the following 
results. 

Theorem 1 Let D be the unit disc inM. 2 . If u : D \ {0} — > N is a W^-map with finite 
energy and satisfies the following equation: 

t(u) = aVu + g 
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where a G C°(D), and g G L P (D,TN) for some p > 2, then u may be extended to a map 
u e W 2 >p(D,N). 

Theorem 2 Let (M,g) be a closed Riemann surface, and N a compact submanifold ofR K . 
Let f be a smooth positive function on M. Assume that u k G W 2,2 (M, N) is a sequence which 
satisfies 

fr(u k ) + V/ • Vu k = a k 

and 

f \Vu k \ 2 fdV g <C, 
Jm 

where a k G L 2 (u k ~ 1 (TN)) and satisfies 

ll a fcllL 2 "~ * 0> as k — ► +00 
If p is a blow-up point of the sequence, i.e. 

limliminf / \Vu k \ 2 fdV g > 0, 

then, p must be a critical point of f . 

Applying the above theorem to the inhomogeneous Landau-Lifshitz equation, we can par- 
tially answer the question mentioned before. Concretely, we come to the conclusions 

Theorem 3 Let (M, g) be a closed Riemann surface, and S 2 a unit sphere with stan- 
dard metric. Suppose that the coupling function f is smooth and positive on M and u G 
L 2 ((0, 00); W 2 ' 2 (M, S 2 )) is the unique weak solution for the initial value problem of the inhomo- 
geneous Landau-Lifshitz equation with initial map uq G W l,2 (M, S 2 ). If u(t) = u(-,t) blows up 
at time infinity, then, the blow-up points must be the critical points of the coupling function f . 

2 Removable singularity 

It is well-known that the removable singularity theorem of Sacks-Uhlenbeck says that a harmonic 
map from D \ {0} —> N with finite energy can be extend to smoothly. The main aim of this 
section is to generalize Sacks-Uhlenbeck's theorem to the present case, i.e., to prove Theorem 
1. The method adopted here is due to Sacks-Uhlenbeck essentially. One still sees that the Hopf 
differential is the key in the proof. However, in our case, the Hopf differential is no longer 
holomorphic, thus the proof should be a little more delicate than theirs. 
Let us first recall the e-regularity discovered by Sacks and Uhlenbeck. 

Lemma 2.1 Suppose that u G W 2 ' 2 (D, N) satisfying 

r(u)=geL 2 (D,TN). 
Then there exits e > such that if J D |Vu| 2 < e we have 

\\u - u\\ W 2,2 {Dl) < C(||Vn|| L 2 (D) + \\g\\ L 2( D) ). 
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Here u is the mean value of u over the unit disc and Di is a disc with radius tt and centered at 

J a 2 

the origin. 

Proof: cf [HUH, or [D],or [EQ- 

Using the standard elliptic estimate, we have 

Corollary 2.2 Suppose that u £ W 2 ' 2 (D,N) satisfies 

r(u)=aVu + g, (2.1) 

where a(x) G C°(D) and g £ L P (D,TN) for some p > 2. Then, there exists e > such that 
whenever f D |Vu| 2 < e we have 

|Vn|(0) < C(\\a\\ c0{D) ,p)(\\Vu\\ L 2 (D) + \\g\\ LP{D) ). 

In this section, we always assume u to be a map from D \ {0} to N which belongs to 
W^(D \ {0},iV) and satisfies the equation (2.1). In order to prove Theorem 1, we need to 
prove the following lemmas. First, we have 

Lemma 2.3 There exists e > such that if f D \Vu\ 2 dx < e, then there holds true 

2— - 

M|Vu|(a;) < C{\\Vu\\ L 2( D2m) + \x\ j> x ||s||z,p(d 2N )) for Vx G Di, 

where C is a positive constant which depends only on e. 

Proof: Fix an xq G Di, we define u = u(x\xq\ + Xq). Then we have 

2 

t(u) = \xo\ 2 g + |xo | VS. 
Notice that |Vu|(0) = |Vit||xo|, then we can get this Lemma from Corollary 2.2. 



Now, let 

^ = {u x ,u x ) - {u y ,u y ) - 2i(u x ,u y ) = 4 ( — , — ), 
where z = x + iy. It is easy to see that 

du du du 

d z y = 8(Au, —) = 8(A{u){du, du) + aVu + g,—) = 8(aVu + 9,q-)- ( 2 - 2 ) 

We need to prove a Stokes type equality for the 1-form z^f. 
Lemma 2.4 There holds true that 

z^dz = / zd^dz A dz. 

z\=r J D r 

Proof: As 

d(z^dz) = dz(z^dz) = zd^dz A dz, 



by applying the Stokes formula, for any rg < r, we have 



/ z^dz = I 

J\z\=r\\z\=m JL 

By (2.2), 



zdz^dz A dz. 

\z\=r\\z\=r J D r \D ro 



i \zd- z ^dz l\dz\ <Cr (\aVu\ 2 + \g\ 2 )dx ^ 



/ z^dz = \^T z 2 ^d6\\ z \ =ro -» 0. 

J\z\=rn Jo 



as tq —* 0. Therefore, to complete the proof of the Lemma, we only need to prove 

z^dz = I 

'\z\=ro 

However the last equality follows from Lemma 2.3. 
Lemma 2.5 There holds 

/ (u r ,u r ) - (ug,ug) = 0(r). 

J D r J D r 

Proof: By a direct computation, we have 

Re(z 2 ^) = -\u e \ 2 + |z| 2 |u r | 2 , 



then 



/Z7T f 
z 2 ^d6\ = \Im j z^dz\ 
J \z\=r 

= \Im zdz^dz A dz\ 

J D r 

< / \zd^dz f\dz\ 

J D r 

<r |aVu + g\ \Vu\dx 

J D r 

<Cr[ (\g\ 2 + \aVu\ 2 )dx. 

J D r 



i.e. 



/ \u r (r,6)\ 2 r 2 d6 - / \u e {r, 9)\ 2 d9 = 0(r). 
Jo Jo 



r /•2ir 



Therefore 

((u r ,u r ) - (ue,ug)) = I I (\u r \ 2 - -^\ug\ z )rd9dr 
D r J o Jo r 

-Oir)dr = / 0(l)dr 
or J o 

= O(r). 

Proof of Theorem 1: As in |S-Uj . we approximate u by the function q which is harmonic 
on every domain, 

D rn (r ) = {z : 2- m -V < \z\ < 2~ m r }, 



and equals to 

r>27T 



W u(2- m r ,9)d9 



and 

r 2n 



1 f /7r 

respectively on the boundaries {z : \z\ = 2~ m ro} and {z : \z\ = 2~ m ~ 1 ro}. Then q is piecewise 
linear in logr and depends only on the radial coordinate. Now, for 2~ m ~ 1 r < r < 2~"V , 

k(r) - u(r, 9)\ < \q(r) - q(2~ m r )\ + \q(2~ m r ) - u(r, 9)\ 

/2tt 
^ |n(2- m r ,^')-^^)l^' 

< C suy 2 -m-i ro <\ r \< 2 -m ro r\Vu\ 
<C(\\Vu\\ L 2 {D2r) +r 2 -l\\g\\ LP{D2r) ). 

Now, we estimate the difference between q and u: 



2tt 

2~ m r 



I |V(n-g)| 2 = J2™ =0 r[\u(r,9)-q(r))(u r (r,9)-q'( 
— / (q — u)A(q — u)dx. 

J D r 

Since q'(r) = constant x £ on 

/ (M(r,0)- 9 (r)) ? '(r)J = O, Vr = 2~ m r , 
Jo 



r))de\t 



Hence, 



r2ir 

TZ=^ I Wr,9) - q(r))(u r (r,e) - q> (r))d9\lZ\ o 



/ZTT 
(u(r ,9) - q(r ))u r (r ,9)d9 
o 

/2tt 
(u(2' m r , 9) - q(2~ m ro))u r (2- m ro, 9))d9. 
o 



By Lemma 2.3, we have 

f27T 



/■ITT 
(u(2- m r ,9) - q(2~ m r ))u r (2~ m r ,9)d9 

< ||n(2- m r , 9) - q(2~ m r )\\ L ^ sup r=2 - mrQ r|V«(r, 9)\ 

-► 



(2.3) 



asm-> +oo. 



Moreover, we have 

/27r / p2n p2n \ 2 

(u(r ,6)-q(r ))u r (r ,6)d6 < r Q I (u(r o ,0) - q(r Q )) 2 dO \u r (r ,6)\ 2 d6 
\J Jo 

< (/ |n e (r ^)| 2 ^J U ^ r 2 K(r ,#)| 2 <ifl 



1 



/2tt 
(|n e (ro,^)| 2 + h r (r ^)| 2 r 2 )^ 




r 

< \\q- u \\l -1 



= ?/ V^o,#)| 2 r o d0. 
^ j 

(2.4) 

and 

Jo — u )(A(q — u))\ = / \q — u\ x \A(u)(du, du) — aSJu — g\dx 

>(.D ro )(PllL°° / |Vn| 2 dx + ^r ||aVu + <7|| L2(Dro) ) 

(2.5) 

Obviously, for any 1 > 5 > 0, we can always pick up ro which is small enough such that 

I \(q - u)(A(q - u))\ <S( I \Vu\ 2 dx + r ). 

Applying Lemma 2.5, we get 

/\V{u — q)\ 2 dx > / {uQ,uo)dx 

= \ \ ((u e ,u e } + (u r ,u r ))dx 
J D rQ 

+ 1 / ((u e , Ug) - (u r , U r ))dx 
J D rQ 

= if \Vu\ 2 dx + O(r ). 

J Dr Q 

Then, from (2.3), (2.4), (2.5) and (2.6) we can derive 

A / |Vn| 2 < r / |Vu(r , e)\ 2 r d9 + Cr , 

JD ro JO 

where A is a positive constant which is smaller than 1. 
Set 

f(r) = [ \Vu\ 2 dx, 

J D r 

then we have 

A/(r) < rf'(r) + Cr. 



(2.6) 



Hence, 



> -Cr-\ 



r 



By integrating the above differential inequality over the interval [r, 1], we obtain 

f(r) < Cr x J s~ x ds + /(l)r A < Cr x . 

By applying Lemma 2.3, it follows from the above inequality that 

\Vu\{x) < |ar| A — 1 . 

Thus, we can complete the proof of the theorem by standard elliptic estimate theory. 

3 A variational formula 

For the inhomogeneous functional Ef(-) defined on W l,2 {M, N), we can see easily that the first 
variational formula at point u G W 2 ' 2 (M, N) can be written as 



dE f (0 = [ < fr(u) + VuVf, Z>dV g , 



for any £ £ T U W 1,2 (M, N). Here, we need to derive another formula of Ef(-) with respect to 
the variation of the domain manifold. The following calculation is essentially due to Price ([P]). 

Take an 1-parameter family of transformations {<f) s } of M, which is generated by the vector 
field X, we have 

E f {uo<f> s ) =~[ \V(uo0 s )\ 2 f(x)dV g 

M 



If ^2\d(uocp s )(e a )\ 2 f(x)dV g (x) 
\ ! \du(^(e a ))\ 2 (Mx))f(x)dV g (x) 
\ ! ^|dn(^,(e a ))| 2 (x)/(0_ s )Jac(^- 1 )^, 



where {e a } is a local orthonormal basis of TM. Noting 

^Jaci^dVg^o = -div(X)dV g , ^/(<M = -df(X), 

we have 

—E f (u o S )| S=O = - 1 / \Vu\ 2 fdiv(X)dV g - ~ / df(X)\Vu\ 2 dV g 
as zj M zj M 

+ E Q / (du(V ea X),du(e a ))fdV g . 
J M 

So, we have proved the formula 

dE f {u)(u*(X)) = -\l \Vu\ 2 fdiv(X)dV g -\[ df(X)\Vu\ 2 dV g 
+ E Q / (du(V ea X),du(e a ))fdV g . 

J M 



$ 



4 The proof of Theorems 



The task of this section is to prove our theorem 2 and 3. In fact, what Theorem 2 concerns is 
just the blow-up analysis for a so-called Palais-Smale sequence of Ef(u). We will focus on what 
occur if the sequence is not compact in the Sobolev space W 1,2 (M, N). 

Proof of Theorem 2: By the assumptions stated in Theorem 2, {u k } C W 2,2 (M, N) is a 
Palais-Smale sequence of maps from M into N. Then, it satisfies 

fr(u k ) + VfVu k = a k (4.1) 

E f (u k ) < C, (4.2) 

where a k G n A T 1 (TA r ) and satisfies 

|K|| L2 ^0. (4.3) 
First, we note that in a local complex coordinates (4.1) can be written as 

fro{u k ) + V / • V u fc = \(3\a k 

where tq and Vo are the operators defined on M? with standard Euclidean metric, since r is 
conformally invariant operator. Without loss of generality, we may assume g = dx 2 + dy 2 on a 
complex coordinate system. 
Set 

S = {x: limliminf / \Vu k \ 2 dV g > 0}. 

Usually, we say x is a bubbling point for the sequence {u k } if and only if x G S. It is easy to 
see that S contains only finite points. By the Lemma 2.1, for any xq G S, we have 



liminf / \Vu k \ 2 dV g > e, for any r > 0. 

fe^+oo J D ( x ) 



r(x) 

By the weak compactness of W 1,2 (M, N) we know that there exists a subsequence of {u k }, 
still denoted by {u k }, and u G W 1 ' 2 (M,N) with 

Ef(u) < +00, 

such that {u k } converges weakly to u in W 1,2 (M, N), which is an /—harmonic map. Moreover, 
Theorem 1, Lemma 2.1, Corollary 2.2 and elliptic estimate theory tell us that u G C°°(M,N) 
and 

u k -> u 

in W x ' q (n, N) for any 17 CC M \ S and q > 1. 

Thus, to prove Theorem 2 we need only to prove 

S C {the critical points of /}. 

Now, pick up a point p G S. As we have pointed out, we may assume g = dx 2 + dy 2 in a complex 
coordinate chart J\f around p. Without the loss of generality, we may assume that p = (0,0), 
Q = [—1,1] x [—1,1] C N and QnS = {p}. Up is not the critical point of /, then, without loss 
of generality, we may suppose that 

df(0) = Xdx, 



9 



where A is a positive constant. Thus, in the neighborhood of p, df(x) = Xdx + 0(r) where 



r 2 = x 2 + y 2 ■ 



We need to choose two functions to cut off the vector field ■§- respectively in x and y 



directions. First, we take a cut off function a G C°° 
where 

5 



which is 1 on [—5, 5], and on [—26, 25] c , 



16 \\ VU \\ 2 C 0(Q)\\J\\C° 

Then, we define the other function as following 

f 1 if \t\ < b' 

{b-t)/(b-b') if b' <t<b 



T)(t) 



(b + t)/(b-b') if -b<t<-b' 
if t > b or t < -b. 



Here b and b' are chosen to satisfy 0<a<b'<b<2a<l where a is a constant such that 



,_ ,o At 

' Vm I < all /II 
2o,2o]x[-l,l] oil " lie 



(4.4) 



Set 



By a direct computation, we have 



d 

X = ri(x)a(y) — . 



div(X) = r]'(x)a(y), 



^2(du k (V ea X),du k (e a )} = rj'(x)a(y) 



dui 



dx 



+ r){x)o-'{y){ 



du k du k , 
dx ' dy 1 



By the formula derived in section 3, we have 

2 



- / r)'(x)a(y)( 
J Q 



du k 


2 


du k 


dx 




dy 



)fdxdy - 2 [ r)(x)a'(y)(-^-, -^)fdxdy 
Jo dx dy 



(A + 0(r))r](x)a(y)\Vu k \ 2 dxdy + 2 / (a k ,u k *(X))dxdy. 

Q 



Note that 



supp(r]'(x)a(y))\Jsupp(r](x)a'(y)) C Q\(-a,a) x (-5,5), 

we can replace Q in the left side of the above equality with Q \ (—a, a) x (—5, 5). 

For arbitrarily fixed a and 5, u k is bounded in W 2 ' 2 (Q \ (— §, f ) x (— |, |)). So, by taking a 
subsequence, we have Vu k — ► Vjj in L 2 (Q \ (—a, a) x (5,5)). Therefore 

2 



/ rj'(x)a(y)( 
JQ 



du k 


2 




dx 




dy 



)fdxdy 



[ V'{x)a{y){ 
JQ 



du 


2 


du 


dx 




dy 



)fdxdy, (4.5) 



and 



)fdxdy 



v(x)a\y){—,—)fdxdy 
-2o,2o]x[-l,l] ^ 



< 211(7'! 

, Ae 

< T' 



C" 11/ NC 



[-2o,2o]x[-l,l] 



|Vu| 2 cfe(iy 



(4.6) 
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where we have used (4.4) in the last inequality. Moreover, once 5 and a are chosen, then 



a,a] X [—S,S] 



\Vu k \ 2 dxdy > e 



when k is sufficiently large. Hence 



(A + Oir))rj(x)o(y)\Vu k \ 2 dxdy + 2 / (a k , u kjf (X))dxdy > -Ae. 



Q 



(4.7) 



In view of (4.5), (4.6) and (4.7), we have 



Q 



Letting b' — > b, we get 



Recall that supper C [— 25, 25] , then 



rj'(x)a(y)(- 


du 


2 


du 


dx 


+ 


dy 


[ 

J |x|=6 


du 


2 


du 


dx 


+ 


dy 



)fdxdy > -Ae. 



)fdy > -Ae. 



5 > 



Ae 



16||Vn||^ 



which contradicts the definition of 5. It means that A must be zero. Therefore p is a critical 
point of /. Thus we complete the proof of Theorem 2. 

Now, we return back to our problem on the location of the bubbling points of the weak 
solutions to the inhomogeneous Landau-Lifshitz equations. Consider the following initial value 
problem: 

f d t u = uA (f(x)Au + V/ • Vu) -iiA(mA {f(x)Au + V/ • Vu)), 
\ u{0) =u {x) e W 1 ' 2 {M,S 2 ). 

Noting \u\ 2 = 1 and the following identity 

u A (u A (f{x)Au + V/ • Vu)) = (u • (f(x)Au + V/ • Vu))u - (u • u)(f{x)Au + V/ • Vu), 
we can see easily that the above equations may be rewritten by 



d t u = f(x)r(u) + V/ • Vu + u A (f(x)r(u) + Vf ■ Vu) 
u(0) = n (x) e PF 1,2 (M, 5 2 ). 



(4i 



Here t(u) = An + |Vn| 2 u is the tension field of map u : M — > S* 2 . 

Tang has ever employed Struwe's method to study the existence and uniqueness of the above 
equation. We outline the argument in [Tj as follows 

1. 3T > 0, s.t. the solution (4.8) solvable on M x [0,T). 

2. u(t) blow-up at finitely many points. 

3. u(t) converges to a u(T) S W 1,2 {M,N) weakly, and on any sub-domain which does not 
contain any bubbling points, u(t) strongly converges to u{T) locally. 

Then, we construct a new flow which stems from ut- Then, by the same argument as in [St] 
we know that There exists T\ > such that the new flow exists on the interval [T, T\ ) and blows 
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up at T\. At each bubbling point, u(t) blows one or more bubbles, i.e one or more non-constant 
harmonic maps. It is well-known that u(t) must lose energy at every bubbling point. Hence, we 
always have a T such that 

f d t u = f{x)T{u) + V/ • Vu + u A (/(x)r(u) + V/ • Vu), 
\ u(0) = u(T) G W" 1,2 (M,S 2 ). 

is solvable on [0, oo). The results in [Tj can be summarized in the following lemma. 

Lemma 4.1 Let (M,g) be a closed Riemann surface and f be a smooth positive function 
on M. For any uq £ W 1,2 {M, S 2 ) there exists a distribution solution u : M x M + — > S 2 of the 
above equation which is smooth on M x M + away from at most finitely many points {xk,tk), 
1 < k < Kq, < tk < oo, which satisfies the energy inequality Ef(u(s)) < Ef(u(t)) for all 
< s < t, and which assumes its initial data continuously in W 1,2 (M, S 2 ). The solution is 
unique in this class. 

It is easy to see the following identity holds for the solution to (4.8) and any < t\ < ti < oo 

E f (u(h)) - E f {u(t 2 )) = - f 2 \\d t u\\ 2 L2 dt. 

Jti 

Hence, it follows 

r+oo 

/ ||^||2 2 <+oo. 

J 

This implies that there exists a sequence dtu(x, U) such that 

\\d t u(x,ti)\\ L 2 -> 0. 

So, {u(ti)} is a Palais-Smale sequence of Ef(u). Therefore, if u(t) does blow up at infinity, then, 
by applying Theorem 2 we obtain the conclusion of Theorem 3. 

As another example, we may also consider the gradient flow of the function Ef, i.e. a solution 

of 

f u t = f(x)r(u) + V/ • Vu, 
\ u(0) = u (x) £ W 1 ' 2 (M,N). 

If u(t) does blow up at infinity, we also have results similar to Theorem 3. 
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